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Abstract 

We study transport in quantum systems consisting of a finite array of 
identical single-channel scatterers. A general expression of the S matrix in 
terms of the individual-element data obtained recently for potential scattering 
is rederived in this wider context. It shows in particular how the band spec- 
trum of the infinite periodic system arises in the limit A — > cxd . We illustrate 
the result on two kinds of examples. The first are serial graphs obtained by 
chaining loops or T-junctions. A detailed discussion is presented for a finite- 
periodic "comb"; we show how the resonance poles can be computed within 
the Krein formula approach. Another example concerns geometric scatterers 
where the individual element consists of a surface with a pair of leads; we show 
that apart of the resonances coming from the decoupled-surface eigenvalues 
such scatterers exhibit the high-energy behavior typical for the 6' interaction 
for the physically interesting couplings. 



Contents 



1 Introductiorj 



2 Serial structure transport 



3 

3 

2.2 Recursive relations for scattering amplitudes! 6 



Preliminaries 



2.3 The S-matrix expression] 7 



2.4 Relation to band spectra of periodic systems] 9 



1 



3 Serial graphs 



10 



3.1 Loop array^ 11 



3.2 Band spectrum of an infinite loop array] 12 



3.3 Comb graphic 13 



3.4 Scattering on a comb graph] 15 



3.5 The resolvent poles] 16 



3.6 Band spectrum of an infinite comb and resonance^ 19 



4 Serial structures of mixed dimensionality 



20 

4.1 Coupling of leads to a surfac^ 21 

4.2 'I'he single-element S-matri^^ 22 



4.3 A "bubble" on the Une 24 



4.4 A "bubble" array] 30 



A A direct derivation of the comb— graph S— matrix 



31 



1 Introduction 



A rapid progress in experimental solid state physics has expanded dramatically the 
list of situations in which consequences of the basic equations of quantum mechanics 
may be tested, since the interaction is prescribed by the experimentalist by the shape 
design of the structure in question, material choice, etc. One of the frequently 
occuring cases in mesoscopic transport involves ge of a quantum particle 

through a serial — or finitely periodic — structure obtained by arraying a certain 
number N of identical scatterers. 

Our aim in the present paper is to study this situation under the assumption 
that the individual scatterers have a single transport mode. For a collection of 
mesoscopic elements connected by quantum wires, this is certainly an idealization. 
We can adopt this approximation provided the transverse modes in the wires are 
well separated and the distances between the scatterers are large enough so that the 
intermode coupling and influence of the evanescent modes can be neglected. 

Such a single-mode transport is often investigated in literature, with the S- 
matrix obtained either "inductively" by adding the scatterers successively or by 
means of the transfer matrix method. It is difficult to collect all relevant references 
but a representative sample is given in ^6]; an extension of this "factorization" 
method to scattering on graphs was proposed recently in |^|. On the other hand, 
the mentioned methods are typically used to evaluate the S-matrix numerically and 
give little insight, say, into its dependence on the number N of the scatterers. To 
this purpose a transparent expression for the S-matrix is needed. 

Such closed-form formulas were derived recently in one-dimensional potential 
scattering, first for an array of 6 interactions |Tl], ]T2], ^ and then for an arbitrary 
finitely periodic potential ||3^ 



45| ; also the number of bound states has been dis- 
cussed in this setting ¥T]. When this work was in its final stage, another analysis 
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of this situation appeared |jTT[ which investigated in detail the distribution of scat- 
tering resonances and the corresponding time delay. Our basic observation is that 
the input for the S-matrix expression are the individual-element scattering data, 
and thus the result can be applied to scattering on an array of arbitrary "black 
boxes". In the next section we rederive the result of [3£, |5| in this more general 
context. 

We apply the whole scheme to two models of quantum particles confined to 
a graph in Sec. 3. The first model consists of an array of planar loops joined 
by a pair of external leads. The system is placed into a homogeneous magnetic 
field perpendicular to the graph plane. From several ways how to couple quantum 
motion on a branching graph in a self-adjoint way we investigate only 6 coupling. 
The next model discussed in this section is made of identical appendices attached 
to a line, called comb-shaped graph. The motion on the line is free, however the 
particle can interact with an external (scalar) potential on appendices. We study 
the most general self-adjoint way of coupling with the wave function continuous 
on the line, the only limitation comes from time-reversal invariance. It leads to a 
three-parametric family of solutions. Dependence of transmission probability and 
resolvent poles on the number of appendices and the coupling parameters is analyzed 
and band spectrum of an infinite comb is discussed as well. 

A step more complicated situation is investigated in section 4; scatterers with 
dimension two are coupled to single-mode leads. First, we impose suitable boundary 
conditions at the junction ensuring self-adjointness and show their relevancy to 
the experimentally observed data. Then we specify our formalism to spheres and 
find explicite expressions for the transmission coefficient. We prove existence of 
resonance peaks and the background which dominates at large energies and decays 
not slower than E^'^ as E oo . Furthermore, we conjecture that the "coarse 
grained" transmission probability (averaged locally over the resonances) has the 
E~'^ decay typical for the S' interaction. Finally, we present band spectrum of 
infinite system. 



2 Serial structure transport 
2.1 Preliminaries 

Consider an equidistant array {Sj : j = 0, . . . , A^— 1 } of identical scatterers placed 
at the line points Xq + ji . The spacing £ > in this convention includes only the 
distance between Sj and not the possible size of the scatterers themselves. 

Let us review briefiy basic notions concerning a single-mode scattering on a 
sole scatterer S placed conventionally at the point Xq = . On the two halfiines 
attached to S the particle moves as free, so the on-shell S-matrix at energy k"^, 
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couples the coefficients of tfie asymptotic solutions 



^b(x) - I + ^-""^'^ ■ ■ ■ < ° (2 2) 



In particular, we have 



t r 
r t 



(2.3) 



where r, t and f, t are the left-to-right and right-to-left reflection and transmission 
amplitudes, respectively. We shall consider only the non-dissipative situation when 
S is unitary; its dependence on the momentum k will be indicated only if necessary. 

Since the interaction responsible for the scattering is localized by assumption, 
solutions to the Schrodinger equation acquire the asymptotic form outside the scat- 
terer. Hence S may be expressed alternatively in terms matrices used in the theory 
of ordinary differential equations. One is the "coefficient" transfer matrix M relat- 
ing the solutions to the left and to the right of the scatterer by 

It is straightforward to see that 

M ^ i ( «:7 ; ) . (2.5) 

and vice versa. 

Another one is the transfer matrix which relates the values and derivatives of the 
two solutions, 

I -'(0+) )-^\ n'{0-) ) ■ (2.7) 
Substituting the boundary values of u{x) = Q<-kx j^^^-ikx f^j, ^ < and u{x) = te'^^ 
for a; > into this relation, we get a pair of equations for r, t which is solved by 

_ L21 + ik{L22-Ln) + k'^Lu 
L21 - ik{L22 + Lii) - k'^Li2 ' 

(2.8) 

^ 2ik det L 

L21 - ik(L22+Lii) - k'^Li2 
In the same way we get the right-to-left amplitudes, 

L21 + ik(Lu-L22) + k'^Li2 

T = 

L21 — ik(L22+Lu) — k'^Lu ' 

(2.9) 

2ik 



L21 - ik{L22 + Lu) - k'^Li2 
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Combining (|2.5| ) with (|2.8|) , ( p.9|) we can express the "coefficient" transfer matrix as 
/ L21 + ik{Ln +L22) - k'^Lu L21 + iA;(Lii -L22) + k'^L^ 



M 



1 
2ik 



\ -L21 + ik{Ln-L22) - k'^Li2 -L21 + ik{Ln + L22) + A;^^^ 
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(2.10) 

The class of admissible transfer matrices is restricted by the S-matrix unitarity. In 
particular, the conservation of probability current, |tp 
I det L I = 1 ; we shall write therefore 



'rP+|t| 



1 implies 



detL 



t 



-: e 



2iip 



(2.11) 



One can also express L by means of the other two matrices. For instance, suppose 
that ( |2.3| ) is given. Using the relation ( |2.7| ) with the boundary values for the left- 
to-right scattering, the same for the right-to-left case with the explicitly written 
and t = t det L , we get a system of four linear equations for Ljk, 



Lii(l+r) + ikLi2{l-r) = t, L2i(l + r) + ikL22(l-r) = ikt 



122(1+ +ikLi2il-f) = t, L2iil+r) +ikLnil-r) = ikt. 



(2.12) 



Only three of them are independent; it is straightforward to see that the first with 
the third, and the second with the fourth equation lead to the same relation which 
is solved by 

(l+r)(l-f) t(r-f) 



^22 



L 
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(l-r)(l+r) (l-r)(l+f) ' 

The same pairs of equations allow us to express L12 and L21 in terms of Lu and 
L22; in combination with the last relation we find 

(1+r) t , ... fl-f) ikt 



ik(l—r) ik{l+r) 



^21 



-ikLii 



ik{l+f) (1 + f) 



We have still the condition (|2.11| ). Computing the determinant with the help of the 
above relations, we get 

t{2Lu-t) _ t 
t' 



detL 



(l-r)(l+r) 

One can express Lu from here and substitute into the formulas for the other ele- 
ments; this yields finally 



2t 



I tt + (l + r)(l-f) 



ik 



- (l+r)(l+f)l \ 



(2.13) 



\ik\tt - {l-r){l-r)\ tt + (l-r)(l+f) j 



The determinant of this matrix is equal to the middle expression of (|2.11|) and 
substituting into ( ^.8] ), ( ^.9[ ) one can check that ( |2.13[ ) indeed represents the inverse 
transformation. 
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Furthermore, the unitarity of S has a stronger consequence. It is well known 
that a general 2x2 unitary matrix can be parametrized by four real numbers as 

gi(a+<s) cos /3 e'^-^-") sin /? \ 
-e*("-^) sin /3 e-*(°+^) cos /? ) ' 



Using this for 5* given by ( ^.31 ) and substituting into ( |2.13| ) we find that 

L = e'^PC , C real with det £ = 1 , (2.14) 
where ip := a + 5 . Notice that M given by (|2.1CI| ) has then the following property: 



Mn = e-2^'^M22 , M12 = e'^^'^Msi . (2.15) 

2.2 Recursive relations for scattering amplitudes 

Before we derive the mentioned closed-form expression, let us recall the usual fac- 
torization technique. We index the transmission and reflection amplitudes for the 
array by iV . In analogy with ( p.3|) we have 

( \ _ I \( 



(2.16) 



where e := e*'^^. Next we add the (A^+l)-th scatterer to the right side of the array 
for which 

[ B^ - [ t e'^f U- ' 



(2.17) 



where, of course, B^ = y4+ and = A^. In analogy with ( |2.5| ) we rewrite the 
last two relations in the "coefficient" transfer matrix form. Multiplying the two 

matrices we get (^^+^ = M(J^jj^ with 



M := 

V i^{-e'''r{tjM-rNrN)-rN) ^(1 - e'^rf^) 

Comparing this with (|2.6| ) we find the sought recursive relations 



-2N~ , ^Ntt ~ Un /0 10^ 

TN+i = e r+ — — , t^+i = . 2.18 

1 — e^^rj^r 1 — e^^rr^- 



If we modify the argument by adding the (A^+l)-th scatterer to the left of the array, 
we get in the same way 



rN+i = e r + , t^+i = =2]v — ' ■ 2.19 

I — e rr^ 1 — e r^r 



Since the "component" S-matrices are unitary, it is straightforward to check by 
induction that the same is true for the total S-matrix. The relations ( p. 181) and ( p. 191) 



have a transparent meaning: expanding the fractions into geometric series we obtain 
expressions containing sums of contributions from various scattering processes. 
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2.3 The S— matrix expression 

On the other hand, the recursive expressions do not relate directly the S-matrices 
of an individual scatterer and that of the whole array. To this end, notice first that 
in view of (|2.15|) we can write the M matrix of the j-th scatterer as 

M, = e^^ [ ^^f^ ^''^ \ , (2.20) 



where e := e as above and 

R := e ^Moo = = — = h i — ; £12 

t t 2 V 2A; 2 , 



(2.21) 



By definition the "coefficient" transfer matrix of the array is obtained by multiplying 
successively the matrices {^^). We denote M^") := M„M„_i . . . Mq . It is easy to 
compute the first few matrices M*^"^; this inspires us to look for the general product 
in the form _ _ 

where the coefficients have to satisfy the recursive relations 

7„+l = C7n + 5n, 5n+l = C^n+ln (2-23) 

with 60 = 1 and 

7«^C:^p, (2.24) 
which follows from A/r("+i) = M„+iM(") . Since 

det M,- = e^^^ - \S\^) = e"'^ ^^t^ = ^^'"^ ' (2-25) 
and consequently, det M^*") = e2*("+i)^|5|2'^+2Q^^|2_ |^^|2^ ^ e2*("+i)'^, we have 



7n = e*''" Vl'^nl' + I'^l"''^-' (2.26) 



with a phase factor to be determined. Substituting into the relations (|2.23| ) we get 



e*^"+^ V^|<5„+i|2+|5|-2"-4 = ^e^e„^|^^|2^|^|^2n-2 +5^^ 

(2.27) 



-2n~2 
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We express e*^" from the second equation and substitute into the first one; this 
yields 

Sn+2 - (C + O'^n+l + (ICr - l)5n = . (2.28) 

Now So = I and 6i = ( + ( , so ( |2.28| ) is solved by 

c + c 



Sr. 



n/2 



where f/„ is the Chebyshev polynomial of the second kind. Since |Cp— 1 = |'S'|~^ 



by ( |2.25| ), its argument can be more compactly written as Re (ei?) . Using ( p.26| 
and (|]23) again, we find M^") = e2*("+i)^A^(") with 



\ 



/ ^n+ig-ie„^i ^ |5|2f/„(Re {eR)y £"5t/„(Re (ei?)) 

^ e"5f/„(Re(ei?)) ^n+ig*e„ ^ |^|2f/„(Re (ei?))2 / 



where 



f/„+i(Re (ei?)) - eRUniRe (eR)) 



(2.29) 
(2.30) 

'l + |5|2[/2(Re(£i?)) 

Now we may employ ( |2.5| ) and ( p.ll| ) to find the sought formulas for the array of 
scatterers; it is sufficient to put n = N — 1 . We arrive at the following conclusion: 

Theorem 2.1 With the given notation, the transmission and reflection amplitudes 
of an N -element serial structure express as 



t 



N 



1 + |5|2[/^_i(Re(ei?))2 
^(^-idN-i SUN-i{R(i{eR)) 
1 + |5|2f/^_i(Re(£i?))2 



(2.31) 



(2.32) 



where the phase factor is given by (^.3(\ ). In the same way the right-to-left ampli- 
tudes are ij^ = t^e-''^^^'^ and 

e^^-'e-''^~^SUN-i{Re{eR)) 

fN — ; • 

^l + \S\mN-i{Re{eR)y 

In particular, the transmission and reflection probabilities are the same in both di- 
rections and equal 

1 , ,o \S\^UN-i{ReieR)y 



t 



N\ 



l + \S\^UN-i{Re{eR))^ 



l + \S\^UN-i{Re{eR))^ 



(2.33) 



Recall that 
\S\' 



Re{eR) = Re 



-i{ke+(p) \ 



t 



1 



Re 



t 



it is obvious from ( |2.33| ) that the probability current is preserved. 
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2.4 Relation to band spectra of periodic systems 

Consider now an infinite periodic array of identical scatterers S joined by line 
segments of length i . The one-period transfer matrix is T = LU^{k) , where 



cos k^ T sin k^ 

k 

-k sin ki cos ki 



coresponds to the segment. The band spectrum of the problem is given by the Bloch 
condition, det (t — e*^^ = , or 

g2ie _ ^ie^^ T + det T = . (2.34) 

In view of ( |2.11|) , detT = detL = e^**^, so the condition may be written as 

e-^rT = 2cos(^-(^) . 

To express the Ihs, we employ ( |2.13|) which yields 

tr T = tr T cos ki + -L^i - kL^ sin ki = — ^ e'''^ + -= e-'^^ , 

\k ) it 

and since 

tt - rf e^*"^ 



t I 
by ( p.5|) and (|2.15| ), we arrive finally at 

Re(£i?) = + — = cos(^-(^). (2.35) 

The //isas a function of k is typically oscillating. Since the amplitude > 1 

unless the single-element scattering is reflectionless, the periodic spectrum has gaps 
in general. 

The relations ( p.33[ ) show how the band spectrum arises in the limit N ^ oo 
of the serial-structure scattering. The Chebyshev polynomials 

,^0 m\{n-2m)\ 

are oscillating within the interval [—1, 1] . The easiest way to see that is to use the 
representation 

Un{x) = — ^T-^-^ , ^ :=arccosx. (2.37) 

Thus Uj\f^i has — 1 roots in [—1,1] and each band contains at least — 
1 points where \tj\f{k)\'^ = 1. Possible additional points with this property can 
come from zeros of the single-element reflection coefficient. Properties of Chebyshev 
polynomials yield also lower bounds to the envelope of the transmission probability 
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oscillations. The representation (|2.37|) implies |f/„(x)| < (1— x^) for |x| < 1, 
and therefore 

^,2 . l-jRejeR)? _ |t(fc)p-(Re(t(A;)e-(^^+^)))2 
'^^ ^' - l + |S|2-(Re(£/2))2 l-(Re(t»e^(^^+^)))2 " ^^'^^^ 

If is a slowly varying function, the rhs reaches its maximum in the 

middle of the band; it is zero at the band edges. However, the transmission can 
vanish within a band only due to a single-element full reflection. This follows from 
another upper bound fll, |t/„(x)| < (n + 1) for \x\ < 1, which yields |tivp > 

On the other hand, to see the behavior of the reflection and transmission am- 
plitudes in the gaps of the periodic spectrum, we have to estimate the Chebyshev 
polynomials outside [—1,1]. By analytical continuation, the relation (p.37|) gives 
Un{x) = sinh ((n+l)arcoshx) / sinh (arcoshx) , so 



n 

Unix) = Y.(x+v^^y 



. n-2k 

ix + Vx^— 1 

fc=0 

Then we have, for instance, the following estimates 

n + x" < Unix) < (n+1) (x + l)" . (2.40) 

The first inequality yields an upper bound, 

\tNik)\^ < o. (2.41) 

' ^ ^' - l + (Ar-l+(Re(£i?))^-i)' ^ ' 

It is clear that |t7v(A;)p = 1 holds only if the same is true for \tik)\'^; in all the other 
cases it behaves as o ((Re(£i?))^^ ^) as A^-^oo. On the other hand, |t^(A;)p = 
holds only if a single scatterer has a full reflection at this energy, otherwise the second 
inequality of (|2.40| ) together with |Re (^-R)| < |t(fc)|~^ and the unitarity relation give 





t(A;)|2^ 


|t(A;)|2^ + iV2 


r(A;)|2(l + 


r(A;)|)2^-2 



\tNik)\' > u.,,.,.^ , . , (2.42) 



3 Serial graphs 

As we have said we want now to illustrate the above results on several examples 
which go beyond the usual one-dimensional potential scattering. In this section we 
shall discuss the situation the serial structure is a graph. 

Although nonrelativistic quantum mechanics for quantum particles confined to 
a graph has been considered already several decades ago in connection with the 
free-electron models of hydrocarbons it became a subject of intense interest 
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\ u, y M 



Figure 1: A loop-graph scatterer in a magnetic field 
only recently as a tool to describe systems of quantum wires — see 0, ^, |10 



III, |T5|, |T6|, [13, [T|, |g, H, |5l 13, m, |3T], g and references therein. There are also 
other systems for which graph description could prove to be useful such as objects 
composed of carbon nanotubes |30| . 

Graph systems are attractive because they are often explicitly solvable; on the 
mathematical level they represent systems of ordinary differential equations in con- 
trast to partial differential equations with nontrivial boundary conditions needed to 
describe quantum wire systems in higher dimensions. At the same time, the sim- 
plified description may still preserve basic features of the real system. Of course, 
replacing a branched waveguide system by its skeleton graph is a nontrivial approx- 
imation, but we shall avoid discussing this point here; some comments can be found 
in 



1^, and references given there. 



3.1 Loop arrays 

In our first example an individual scatterer is a planar loop C with a pair of external 
leads placed into a homogeneous magnetic field — cf. Fig. 1. We suppose that the 
field is perpendicular to the graph plane. The corresponding vector potential in the 
circular gauge is A{x) = ( — ^By, ^Bx, Oj and 



c Li+ L2 Li+ L2 

where 5* is the loop area, Lj are the lengths of its two branches, and $ is the 
corresponding magnetic flux. 



Such a scattering system has been considered earlier — see p, 0, ^ and refer- 
ences therein. The Hilbert space for the graph of Fig. 1 is the orthogonal sum of four 
spaces referring to the graph hnks; its elements will be denoted as (/i, mi, U2, /2) 
with the coordinates at the loop taken anticlockwise. We suppose that the particle 
has a unit charge and apart of the magnetic field it moves as free on the graph. There 
are various ways how to couple the operators {—idx — A)^ from different graph links 



in a self-adjoint way |2g, |3J]. For the sake of simplicity we restrict ourselves to 
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the usual S coupling, i.e., we impose the boundary conditions 

/i(0) = Mi(0) = U2{L2) -f[{0) + u[{0) - u',iL2) = ai/i(0) , 

/2(0) = U2{0) = u,{L,) /^(O) - u[{L,) + ^'2(0) = aMO) , (3.1) 

where the function and derivative values mean the appropriate one-sided limits. On 
the other hand, we do not suppose in general that Oi, 02 and Li, L2 are the same. 

Substituting into ( p.l|) the boundary values of the generalized eigenfunctions 
on the loop, Uj{x) = Aje~'^^^~''^^ + Bje'^^^^'^^^ , j = 1,2, we get a system of six 
equations. Eliminating from here Aj, Bj we arrive at the relations 

/2(0) = ^{{za^B + 2kC)m+tBf[{0)) 

/2(0) = ^({ta,a2B + 2ia^ + a2)kC -Aik''E^E+)fi{0) + {2kC + ia2B)f[{0) 
IkD ^ 

determining the transfer matrix, where 

B = {el-^- — el-){e2-—e2+) = 4e~** sin fcLi sin fcL2 , 
C = -61^62- + ei+e2+ = 2ie"** sin/c(Li + L2) , 
D = £2+ - £^2- + £2+£2-(^i+-£:i-) , 

^'S>/2^±^kiL,+L,)/2. ( ±k{Lr + L2) - ^ \ 



E± = I — ei±e2± 
e^-i = e^(-^^^)^^- , J = 1,2. 
The reflection and transmission amplitudes are given by 

-iaia2B + k{{a2-ai)B - 2{ai + a2)C) + ik'^{4:E_E+- B) 



r{k) 



m 



iaia2B + k{a2+ai){2C + B) - ik'^{AE_E++ B + AC) 
Aik'^De-^''^ 

iaia2B + k{a2+ai){2C + B) - ik^{AE^E++ B + AC) ' 



(3.2) 



respectively. Since B, C, E± as well as are 2'7r-periodic functions of the magnetic 
flux, the same is valid for the reflection and transmission probabilities. Recall that 
if we put e = h = c = 1, then 27t is the magnetic flux quantum in these units. 



3.2 Band spectrum of an infinite loop array 

We illustrate the relation of transmission probabilities of a finite array of loops and 
the spectrum of the corresponding infinite system on Fig. 2. As already mentioned 
|tiv(A;)p is a 27r-periodic function and because the condition (p.35|) determining the 
band spectrum of the infinite system can be reformulated as < |Re(ei?)| < 1 
we show dependence on $ only in the range [0, 27r] (|Re(ei?)| is also a 27r-periodic 
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function). We choose loops with different Li,L2 and ai,a2- Even for a relatively 
small number of loops N = 6 the values of \tN{k)\^ are clearly nonzero in areas 
of parameters $ and k where there are bands of the infinite system and negligible 
where there are gaps. 



3.3 Comb graphs 

Our next example concerns the case of a comb-shaped graph, i.e., a line with a finite 
number of identical appendices attached to it at equally spaced points. Such 
systems have been discussed recently |^2[ following earlier studies of a single-stub 
waveguide 0, 0. 

Comparing to the previous work and the preceding example, we shall discuss 
comb-shaped graphs in a more thorough way. First of all, instead of the 5-coupling 
used above (or the Griffith's boundary conditions in the terminology of 



we 



allow for the most general self-adjoint way in which the stubs can be attached to 
continuous wavefunctions on the line. This amounts to imposing at the junctions 
the boundary conditions adopted from an earlier treatment of the T-shaped graph 
Should the Hamiltonian be time-reversal invariant, the junction is then charac- 



terized by three real parameters. In this framework we are able to handle imperfect 
contacts |]T^; moreover, it is straightforward to modify the results derived below to 
graphs with a ^'-coupling which corresponds to the situation where the junction it- 
self represents a complicated geometric scatterer (see [0, ITB] , more about that will be 



said in the next section). Computing the S-matrix we also assume that the particle 
is under influence of a potential on the stubs; this makes it possible to investigate 
how the band-form zones of high transmission which arise for ^ 1 change when 
an external field is applied. 

The formula derived in the previous section allows us to express the transmission 
and refiection probabilities. In addition to them, one is able to generalize the result of 
23 to the present situation and to find an explicit Krein-formula expression for the 



resolvent of the comb-graph Hamiltonian. This allows us to study their resonance 
structure of the problem which arises from perturbation of the disconnected-stub 
discrete spectrum embedded into the continuum of the line motion. The mentioned 
expression yields an equation from which the resolvent singularities on the second 
sheet can be found. 

After this introduction, let us describe the model. For a greater generality we 
suppose first that the appendices are not necessarily identical. The graph will 
therefore consist of a line with a finite sequence {(s — l)£}^i of points at which 
appendices of finite lengths are attached (see Fig. 3). The state Hilbert space 
of the problem is then H = L^{Tn) := L'^{1R) © (©^^ L2(0, L,)) ; we shall write 
its elements as columns 

( ^ ^ 

Ml 

^ = . 

\un J 
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L 



N 



Figure 3: A comb-shaped graph 



We suppose that the motion at the backbone hne is free while the particle is exposed 
to potentials Vs on the "teeth" ; hence the Hamiltonian of a nonrelativistic particle 
of mass m = 1/2 living on r^v acts as 

{HijUx):=-f"{x), (//^),+i(x):=(-<+Kws)(x), s = l,...,N, (3.3) 

out of the junctions. To make it a self-adjoint operator one has to choose again 
properly the boundary conditions which couple the wavefunctions at the branching 
points of the graph in such a way that the probability current conservation at the 
vertices is conserved. As we have said, we shall require continuity of the line com- 
ponent f of ip at the junctions; to keep a manageable number of parameters we 
assume the Dirichlet conditions at the end of the stubs. Consequently, the domain 
of the Hamiltonian consists of all ip e H with / G AC'^{IR) and Ug G AC'^{0,Ls) 



satisfying the conditions 



/((,_!)£+) = /((.-1)£-) =: /((.-IK), 
UsiO) = 6J((s-l)£) + cX(0), 
/'((,_!)£+) _/'((,_!)£_) = dJiis-m-bsu'M, (3.4) 

Us{Ls) = 

for s = 1, . . . , . In general, bg may be complex, however, we restrict from the start 
to Hamiltonians which are time-reversal invariant and suppose that the coefficient 

matrices Kg = (^J^ '"^ ^ ^^e real. The s-th appendix is decoupled from the line 

by putting 6^ = ; it is then described by the operator h^^ '■~ " + specified 
by the decoupled condition 

Ug{Q) - c,<(0) = (3.5) 

at the junction. 

In what follows we concentrate on the finite periodic case where Lg = L , Vg = 
V , and IKg = IK ; the operator specified by the boundary conditions ( p.4| ) will be 
denoted as = Hn{]K, V) . The potential V is supposed to belong to L^{0, L) . 
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3.4 Scattering on a comb graph 

To write the scattering matrix we need some notation. Let be the unique solution 
to the appendix Schrodinger equation, 

- ul + Vul = k'^UL (3.6) 

with the energy k"^ and the normahzed Dirichlet condition at the outer endpoint, 
ul{L) = 1 — m'^(L) = . If all the junctions are described by the parameters b, c, d 
in (|3.4|) , we put 

/3 := 7^ + ^ (^^] (0), C := cos ki + (3 sin ki, (3.7) 
2k 2k \cui^ — ul) 

and denote again e := e^^^. For a single stub, s = N = 1 , we insert a general- 
ized eigenvector of the form into the boundary conditions ( |3.4D and express 
/(0+), /'(0+) by means of the left-sided limits; this yields the single-element trans- 
fer matrix which can be with the help of ( |3.7| ) written as 




so R = l + ij3 and S = i(3 . Furthermore, ReieR) = (. Using then the explicit 
expression ( p.30| ) of the phase factor together with (|2.31|) and ( p.32| ), we may write 



the left-to-right scattering amphtudes as 

" eUN-2iC) - {W^)Un-i{C) ' eUM-2{0 - (1+^/5) t/^-i(C) ' ^ ' ^ 

In particular, the reflection and transmission probabilities are 

|t^(A;)p = 1 - \r^{k)\^ = -— -i . (3.9) 

1 + P^Un-i{Q^ 



In order to illustrate the usefulness of the formulas ( p.31|) -( ^.33| ), we present in 
Appendix A two other ways to derive the relations 



Let us now illustrate how the the transmission probability depends on the num- 
ber of the teeth and the parameters of the junctions. Fig. 4 shows a typical example 
of the situation. For a certain set of parameters b, c, and d we compare |t(A;) p for one 
and seven appendices with an infinite array of them. The band spectrum is marked 
by thick lines above. The suppression of transmission coefficient at those values of k 
where the infinite system shows gaps is evident even for a relatively small number of 
elementary scatterers, N = 7. In each band we distinguish — 1, i.e. six values of 
ki when |t(A;j)p = 1 as it should be - cf. (|2.37| ). In addition the first and third band 
contain an additional value of ka for which |t(/ca)p = 1 and these values coincide 
with the values of k for which a single appendix allows a perfect transmission. The 
values of k for which the refiection is total fall into gaps of the infinite system as it 
should be {cf. the end of sec. 2.4). 



15 



3.5 The resolvent poles 



The energy dependence of the transmission probabihty through the comb graph is 
related to its resonance structure which in turn comes from perturbation of embed- 
ded eigenvalues of Hn{0, V) referring to the stubs by the coupling between the line 
and the "teeth" . The resonances are conventionally associated with the poles of the 
analytically continued resolvent; we are going to show now how they can be found. 

As in the particular case = 1 discussed in one can derive an expression 
for the resolvent {Hn—z)~^ of the operator Hn{1K, V) by means of Krein's formula 
Appendix A]: it is an integral operator with the kernel 



2N 

E 

m,n=l 



(3.10) 



in this relation Hq is a suitable comparison operator and the vectors satisfying 
Fm{x, z) = Fm{z) belong to the the deficiency subspace (for a given complex energy 
z ) of the maximum common restriction of Hn and Hq . 

As usual we put k := ^/z . We choose for Hq the fully decoupled operator with 
bs = Cs = ds = , s = 1, . . . , N . Its resolvent can be written down explicitly; its 
kernel is a (A^+1) x (A^+1) matrix 



( Riix,y;z) 






M^(«0,l,«L,l) 








V 







(3.11) 

where Vr(tio,n, Wi,n) is the Wronskian of the solutions Mq,™ and to — m" + 

VnU = k'^u specified by the boundary conditions Mo,n(0) = and UL,n{Ln) = 0, 
respectively. The symbols a;< and Xy mean respectively the smaller and larger of 
the variables x, y , and 



Ri{x,y;z) :-- 



Jk\x-y\ 



2k 



(3.12) 



is the free resolvent kernel on the line. A natural choice of the deficiency vectors is 
the following: 



Fnix) :-- 



( Ri{x, {n-l)i;z) \ 




/ \ 

^iL,n(x, Z) 



\ 







(3.13) 



for n = 1,...,A^. To find the coefficients Xmn "we use the fact that for z of 
the resolvent set the resolvent maps L'^^Tn) onto the domain of . Hence the 
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components of the vector at the Ihsof 



( f \ 

Ml 



( 9 \ 



(3.14) 



have to satisfy the boundary conditions (3^) for arbitrary (square integrable) func- 
tions g,vi, . . . ,vn ■ For the sake of brevity we introduce 



K-= Ri{y,{n-l)L;z)g{y)dy, K+n 
Jm 



Un{y,z)vn{y) dy 



with n = 1,...,N; they represent 2N independent quantities. By direct compu- 
tation we find then the needed boundary values from p.lU| )-( |3.14] ), 



<(0) 



2N / N 



dRi{x,y;z) 



M dx 



9{y) dy + 



2N 



s-2 



N 



+2 E^H -E^n^^'""+ E 

^ * 71=1 n=s—l 



1=1 

dRi{x,y; z) 



M dx 



9{y) dy + 



x=(s-l) 



2N I s-1 N \ 

l=\ \ n=l n=s / 



2N 
k=\ 



1 



27V 



(3.15) 



where s = 1, . . . ,N and e := e^^^ . Substituting from here to (|3.4| ) and using the 
fact that hn are independent we arrive at the following system of equations 



N 



ib 



E + (c.<.(0) - «L,s(0)) A;v+M 

n=l '^'^ 



n=l 



-6s I = S 

l = N+s 
otherwise 



(3.16) 



/ = s 
I = N+s 

ML,s(0) 

otherwise 



-4 

bs 



17 



with s = 1, . . . ,N, I = 1,..., 2N . If we arrange the sought coefficients in the 
following way, 



•^1,1; • • • ) '^2Af,l; -^1,2; • • • , A 



27V, 2; 



^1.2Nj • • • , a 



2N,2N 



the matrix of the system (|3.16|) is block diagonal with 2N same 2N x 2N blocks. 
This is important because we are looking for singularities of the resolvent ( p.lOD 
which occur at the points where the coefficients Xmn are singular. The latter are 
solutions to the linear system (|3.16|) , and therefore they have the same denominator; 
looking for its zeros it is sufficient to inspect the determinant of a single block, 



ibs g-.\s—n\ 
2k ^ 



( 



CsU'^s -Ul,s) (0) 5s,N-n 



For a non-real z we have 

[csu'l^, - UL,s) (0) ^ and &sm'l,.(0) ^ 

for all s = 1, . . . ,N , since otherwise the decoupled stub Hamiltonian with the 
appropriate boundary condition at x = would have a complex eigenvalue which 
is clearly impossible. Vanishing of the above determinant is then equivalent to 



Dn := det (z/^^el^""! + 6s,n) = 



where 



dn bl 



u 



L,n 



■- ZIl ZIL I 

2k 2k \cu'L,^-UL,r. 



(0) 



Up to now the stubs could be mutually different. We shall evaluate the determinant 
in the finitely periodic case, I3n = (3 for ra = 1, . . . , . Then can be written as 



D 



N 



[i(3{l-e^) + l)bN-i- 



where 



= i(3 + l, Di = l + 2ip- p^l-e^) , 
and the quantities satisfy the recursive relation 

eDn - 2zDn~i + ebN-2 = . 
Comparing this to the relations between Chebyshev polynomials we find 

Dn = -£^-M^f^iV-2(C)-(l+^/5)t^7V-l(C)) , 



(3.17) 
(3.18) 
(3.19) 

(3.20) 



since ( |3.17| ) is a linear difference equation of the second order with constant coeffi- 
cients which has for the given initial conditions (|3.18| ) a unique solution. 
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We see that DN{k) = if and only if the same is true for the denominator of the 
reflection and transmission amplitude. Hence the poles of the continued resolvent 
and of the S-matrix generically coincide. It might be that some of them vanish 
due to a zero in the numerator, however, this can happen in isolated points only 
because the numerators are analytical functions of the coupling parameters, and 
there is a one-to-one correspondence between the pole trajectories with respect to 
the coupling constant and the embedded eigenvalues of the decoupled Hamiltonian 
H{0,V). 

3.6 Band spectrum of an infinite comb and resonances 

We discuss only the case of free appendices here, i.e. Vg = 0. Investigation of non- 
zero potentials represents no complication in general conception. Typical spectra 
of an infinite comb are presented on Fig. 5. The parameters of such a structure 
could be divided into two groups, viz. i, L and b, c, d. The dependence on i is rather 
simple; a change of i results in scaling of spectrum in k. This is the consequence 
of the fact that the transfer matrix T depends only on the product ki. This is the 
reason why we show all the spectra for i = 1 only. The dependence on L is more 
complicated as can be inferred from Fig. 5. 

A striking feature of such a spectrum are sudden transitions of bands to gaps or 
vice versa at k£ = nvr, n G Z for almost all L. This can be understood, if we write 
down the Bloch condition explicitly as 

f3{k) sm{k£) + cos(H) = cos(^). 

Clearly, this is always fulfilled for ki = nn, unless (3 has a singularity at nn/i. 
Therefore k = mr/i belongs always to the band in this case. The situations when 
j3{mi /€) is infinite requires a further analysis. An effect of a small change e'mk can 
be estimated from I3{k + e)e£ + 1 = (— 1)" cos(6'). With exception of points where 
(3{n7T/i) = 0, small change of e does not change the sign of (3, while the change of 
sign of e leads to the change of sign of P{k + e)e£. Only at vicinity of points where 
P{mT/i) = can the sign of (3{k + e)e£ rest unchanged, cf. Fig 5. 

Of the three parameters determining the coupling of an appendix to the back- 
bone line i.e. b,c,d, we find b as the most suitable to begin with. Putting 6 = 
switches effectively off the appendices as already mentioned in sec. 3.2. This implies 
that we have eigenvalues /c^ embedded in the continuum, where fc„ are solutions of 
(|3.5| ). These values of momentum play an important role from three points of view. 
First, knS are the values for which ml(0) = and the number of zeros of u increases 
by one when k passes to higher values; this means that the whole IR^ is divided 
into disjoint intervals by these values (for a given L). On the other hand, a band 
can belong to two intervals (although this happens only for some L) and states of 
the same band can be described by wave functions with different number of zeros 
of m's. Also, and this situation is far more frequent, one interval can contain two or 
more band, so that the number of zeros of u is not directly related to the ordinal 
number of a band. 
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In order to reveal the second role of these values knS, we have to return to finite 
number of appendices attached to the line. If we begin with one appendix only, 
direct calculation shows that t{kn) = and r(fc„) = —1. Therefore, regardless of the 
number of appendices the full refiection takes place, i.e. tN = 0,rN = —1. 

We close this section with discussion of scattering resonances of our system. 



Recently, it was shown |]TT| that there is certain "band" structure in the spectrum 
of resonances and that this spectrum converges to the energy band spectrum of the 
infinite periodic system in the limit of an infinite number of scatterers. The authors 
also show that one should expect N—1 resonances in each band if the system consists 
of identical cells. These conclusions are in accordance with our results but for one 
difference - the number of resonances in a band is 2N — 1. It is a consequence of the 
fact that the origin of resonances is twofold here. N — 1 resonances comes from the 



spacial order of our periodical structure as in |]1T| and the remaining ones have 
their origin in /c^, which is an A^-fold eigenvalue for the case of 6 = 0. As |6| grows 
this degenercy is lifted and we find A^ resonances in the vicinity of k"^. This is the 
third role of The A^ — 1 resonances related to the spatial setting travel from — oo 
and for small b these two sets are well separated as is shown on Fig. 6. 

There are situations when a resonance or a bound state appears at k = 0. The 
manifestation of this instance is that |t(0)p = 1. This can take place even in the 



case of one appendix. The expression (5) of with explicit form of ul(0) shows 



that this can happen as soon as the condition 

d{c + L) + 6^ = 

is fulfilled. Provided that this condition holds, there is perfect transmission at A; = 
for any number of appendices, z.e.tAr(O) = 1. Besides this there could be other 
situations when tAr(O) = 1 for A^ > 1. If there exists such an n = 1, . . . , A^ — 1 that 

/nvrX ^ (dc + dL + b^)i 
llv j = ' + 2(c + L) • 

we have again 1^(0) = 1. This corresponds to Un~i{C{^)) = in (|3.9| ). Direct 
inspection confirms that the wave function does not belong to L^^Tn) and that 
these states are resonances (and not bound states). 



4 Serial structures of mixed dimensionality 

In this section we want to treat the situation when the scatterers connected by 
single-mode leads have a higher dimension. For the sake of simplicity, we shall 
consider only the simplest possibility when the dimension is two, i.e., the scatterer 
is a surface. Such systems can be realized in both the sohd state (recall, e.g., the 
"bamboo defects" in nanotubes [^) and electromagnetism (for fiat resonators), 
but we avoid discussing examples and the conditions under which these models are 
realistic. We suppose that the surface is smooth, bounded, and connected, with or 
without the boundary. Although it makes no difficulty to let the particle on the 
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surface interact with an external potential field, we will regard it as free, i.e., its 
Hamiltonian will be (in appropriate units) just the corresponding Laplace-Beltrami 
operator. 



4.1 Coupling of leads to a surface 

The basic question for the described serial structures is the way in which the leads 
are coupled to the scatterers. The physical condition is again a conservation of the 
probability current, which translates into the self-adjointness requirement of the 
corresponding Hamiltonian. Since the coupling is local, we may disregard geometri- 
cal peculiarities of the lead and the surface and consider the setting when a halfiine 
is attached to a plane. The state Hilbert space is then L2(iR-) © L'^{IR^) and the 
Hamiltonian acts on its elements (^^^^ as make it self-adjoint one has 

to impose suitable boundary conditions which couple the wavefunctions at the junc- 



tion. A general solution to this problem is given in Ref. |21|. The conditions are of 
the form 

<l)[{0-) = AMO-) + BLo{(j)2) , Li{<f)2) = CMO-) + DLoih) (4.1) 
together with several "exceptional" classes, where 

Lo(02) = hm , 

r^o+ inr 

Li(02) = lim (02(^) - -^^0(^2) Inr) 

r— >0+ 

with r := |x| are the generalized boundary values in the plane, and the coefficients 
A, B, C, D depend on four real parameters; evaluating the boundary form, it is 
straightforward to see that they satisfy restrictions 

A,DeR B = 27iC. (4.2) 

A disadvantage of this result is that it tells us nothing about physical relevance of 
the coefficients values in the boundary conditions ( [4.1[ ). The choice of the coupling 
depends on particular properties of the junction it models, of course, but one would 
like to select a subclass representing a "natural" coupling. One way to achieve this 



goal was suggested in [23]: comparing the scattering matrix of the junction given by 
(|4l| ) with the low -energy behavior of scattering in the system of a plane to which a 
cylindrical "tube" is attached, and taking into account the condition ( |4.2|) , we arrive 
at the identification 

where p is the contact radius. Physical relevance of these conditions was illus- 
trated in |2^ by explaining the experimentally observed distribution of resonances 



in a microwave resonator with a thin antenna. Motivated by this, we will use in 
the following ( [4.1[ ) and ( [4. 3] ) to describe the coupling between the leads and the 
scatterers. 
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4.2 The single— element S— matrix 

Using the local character of the boundary conditions derived above we apply them 
to coupling of a pair of halfline leads to an arbitrary surface G . The only restriction 
is that that the junction may not belong to the boundary of G if it has any. We 
shall compute the transfer and scattering matrices for such a system. 

As we have said the Hamiltonian is a Laplace-Beltrami operator on the state 
Hilbert space L'^{G) of the scatterer. We shall characterize it by its Green's function 
G{., .; k) , i.e., the integral kernel of its resolvent which exists whenever fc^ does not 
belong to the spectrum. Its actual form depends on the geometry of G but we 
shall not need it. What is important is the character of its singularity. As a smooth 
manifold, G admits in the vicinity of any point a local Cartesian chart and the 
Green's function behaves as that of Laplacian in the plane, 

G{x,y;k) = - — \n\x-y\ + 0{1) , \x-y\^0. (4.4) 
2n 

Looking for scattering solutions to the Schrodinger equation, we need a general 
solution to the Laplace-Beltrami equation on G for the energy k'^. Without loss of 
generality, we may write it as 

u{x) = aiG{x,xi; k) + a2G{x,X2; k) , (4.5) 

where xi,X2 are two different points of G at which the leads are attached. The 
generalized boundary values (labelled by the point at which they are taken) of this 
solution are then 

Lo[xj] = -TT ^ Li[xj\ = aj^{xj,k) + a3_jG{xi,X2;k) (4.6) 

ZTT 

for j = 1,2, where 

In \x—x 



e(x,; k) = \im 



G(^Xj "X j ^ ^) ~\~ 



27r 



(4.7) 



Next we denote the wavefunction on the j-th lead as Uj . For simplicity we 
use the abbreviations Uj, u'j for its boundary values; then the boundary conditions 
§3) yield 

/ A Biai Diai 
= Aiui — , ai4i + 025- = CiMi — , 

ZTT /TT 

U2 = -A2U2 + — , a2^2 + aig = C2M2 — , 

where g := G{xi, X2]k). In the the first equation of the second pair we have changed 
sign, because the second lead is identified with M'^. It is straightforward to rewrite 
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these equations as a linear system with the unknown U2, u^-, cti, ^2 and to solve it; 
this gives in particular the transfer matrix, 



9C2 



A1A2 



-A2 



B2C2Z1 

Bi 



where Zj := ^ + and A := g'^ — Z1Z2 ■ Using ( |4.2| ) we find easily 

, , r B2C1 C2C1 

deti^ = -- 7^ = -^TT. (4-9) 

SO det L = 1 if the junctions are identical or the coefficients Cj are real. The second 
possibility if the couplings are invariant with respect to the time refiection, which 
we shall suppose in the following. In that case the transfer matrix simplifies to the 
form 



^ A 



-2 



in particular 



C2 



A{Z^ + Z2)-§A Aa + Zi 



Zi + ttA 



-27rpA \ 



(4.10) 



(4.11) 



^\i{l^-Zi-Z2-^A) Zi + vrAy 

for the physically most interesting class of couplings ([4.3|) . The S-matrix of our 
geometric scatterer is then given by the relations (|2.31|) and (|2.32|) ; in the case 
(|4.11|) we have 

ttA + Zi + Z2 - TT-i + 2ikp{Z2-Zi) + 47r/tVA 



r{k) = — 



vrA + Zi+Za 



2iA;p(Zi + Z2 + 27rA) -47rA;VA 
Aikpg 



(4.12) 



ttA + Zi+ Z2- TT-i + 2iA;p(Zi + Z2 + 27rA) - Ank'^p'^A 



To make use of these formulas, we need to know Zi, Z2, A as functions of 
the momentum k . By assumption the manifold G is compact, so the spectrum 
{^n}'^=i of the Hamitonian is purely discrete and the corresponding eigenf unctions 
{0(x)n}^i form an orthonormal basis in L^{G) . The usual Green's function ex- 
pression then gives 



9ik) = E 



n=l 



(Xi)0„(x2) 



(4.13) 



To express the remaining three values we have to compute the regularized limit 
(|4.71 ). Expanding the logarithm into the Taylor series, we can rewrite the sublimit 
expression as 



G{xj + y/en, Xj] k) + 



Iny/g 
27r 



E 

n=l 



>„(Xj + ^n)(t)n{Xj] 

K - k^ 



Aim 
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where n is a unit vector in the local chart around the point xj . Unfortunately, 
interchanging the limit with the sum is not without risk since the latter does not 
converge uniformly. To see that the result may indeed depend on the regularization 
procedure, it is sufficient to replace y/e by Cy/e at the Ihs. To make idea about this 
non-uniqueness, let us compute the difference 



C{xj,k) -C{xj,k') = lim ^ 



(j)n{Xj + y/en)(l)n{Xj) 0„(Xj + ^/en)(f)n{Xj 



This sum is already uniformly convergent, because by standard semi classical esti- 



mates |3^, XIII. 16] the sequence {||0n||oo}^i is bounded with our assumptions and 



Xn = 4:7i\G\ + 0{1) as n — > oo , so 

1 11 



2 ' 



A„— k"^ Xn — k"^ n 

and therefore 



(4.14) 



From the same reason 



makes sense and ^{xj, k) — ^(xj, k) is independent of k . We have therefore 

«..^)^|(^-ji^).c(G). (4.16) 

The constant depends only on the manifold G we will neglect it in the following, 
because its nonzero value means just a coupling constant renormalization: Dj has to 



be changed to Dj+ 2ttc{G) . For a fiat rectangular G we found in an agreement 
with the experiment using c{G) = . 

4.3 A "bubble" on the line 

To make the above consideration more concrete, we shall concentrate in the rest 
of this section on a single example. We are going to consider the case when G 
is a sphere of radius R with the leads attached at the poles, which is the system 
proposed by Kiselev |3^. The most important result of this paper was that apart of 
the resonances coming from the bound states on the sphere, such a scatterer has the 
high-energy behavior similar to that of the so-called 6' interaction |^, |T^, i.e., 



the transmission probability decays as E ^ for E ^ oo . This was established in 
Ref. up to a logarithmic correction. The difference here is that we shall use the 



physically interesting coupling ( |4.3| ), while |Q employed a two-dimensional subset 
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of the conditions (|4.1| ) disjoint with the above one. This leads to a different S-matrix, 
and while the indicated high-energy behavior remains preserved, the argument used 
in to demonstrate it has to be changed in numerous places; this is why we present 
its modified version here. 

The sphere Hamiltonian is chosen in the standard way. Using the spherical 
coordinates, we write it as 

- cot_^^ 1 

with the usual domain. For the sake of simplicity we put i? = 1 in the following; the 
results for a general R can be obtained by a scaling transformation. The spectrum 
of Hg then consists of the eigenvalues are A/,m = ^(^ + 1), / = 0, 1, . . . , m = 
—/,—/ + !,...,/ of multiphcity 2/ + 1 to which the eigenfunctions 



(2Z + l)(/-|m|)! 
^ 47r(/+|m|)! 



P,''"'(cos^)e*"^ 



correspond. The junctions Xi, X2 we place at the points 9 = and 9 = n , 
respectively, where 

so only states with m = can be coupled to the leads. To express the S-matrix, 
we need the quantities 

ZA^) - Z(k) -.^It - E p^) - If + . (4,19) 

where p is the junction diameter which is supposed to be the same for j = 1,2. 
The relations (|4.12|) yield, in particular, the transmission probability in the form 

tik) = ^mkpgjk) 

^' l-27TZ{k){l + 2tkp) -'K^A{k){l + 2tkpy' ^ ' 

where A{k) = {g{k) - Z{k)){g{k) + Z{k)) . 

To prove that (|4.2(]| ) behaves as indicated above at large values of E = k"^, we 
need several auxiliary results. 

Lemma 4.1 The function Z{-) + {—iyg{-) is strictly increasing on the intervals 
(/(/-I), /(/ + !)). 

Proof: We have 

z'{k)H-i)'9\k) = 2kY.^-^-J±-^^ > 0. I 
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Lemma 4.2 We have 
9{k) = ^ 



2/-1 



2/ + 1 



47r 



/(/-I) + 



/or E G (/(/—I), /(/+!)) anc? any positive integer I with the error term independent 
of I . Moreover, there is a K > such that g{k) > K holds for all k large enough. 



Proof: The error can be estimated explicitly; we shall show that 



9ik) - 



-I] 



21-1 



47r 



+ 



21 + 1 



l{l-l)-E l{l + l 



E 



< 



27r 



(4.21) 



To this end we have to find a bound to ( [4.18| ) on the interval (/(/ — I), Z(Z+1)) with 
the two singular terms removed. The terms to the left and to the right of this pair 
form alternative sequences with the decreasing modulus, and as such each of them 
may be estimated by the (modulus of) the first term of such a sequence, i.e., by 



1 

47r 



2/-3 



{l-l){l-2)-E 



and — 

47r 



2/ + 3 



{l + l){l + 2)-E 



respectively. Taking the maxima of these expressions and summing them we arrive 
at ([4.21| ). To get the second claim, we have have to compare this result with a 
lower bound to the modulus of the sum of the two singular terms. The minimum 
7r~^+ 0{l~^) of the latter is reached at E = l"^ — ^ + 0{l~'^) , so this part wins over 
the other one once the error terms become small enough. | 

Lemma 4.3 For a positive integer I and E G (/(/ — !),/(/+ 1)) we have 



1 



21-1 



+ 



21+1 



A'K\l{l-l)-E l{l + l)-E 
where the error term is independent of I . 



(4.22) 



Proof: We again split the two singular terms in ( ^4.22| ) and write the rest as Z_{k, l) + 
Z+{k,l) - +c(G), where 



1-2 



Z_{k,l) 

Z+{k,l) 
It is easy to see that 



E 

•=i 

oo 

E 

j=i+i 



2j + l 



2j ^ 

-Y — — 

47r \ + -E ^0 



1 

A-n 

F = — 



23 + 1 



2i 

-E 



j{j + l)-E ^^f+n+l 



26 



as J — ^ cxD , so the j-th term in Z±{k, I) can be estimated from above by 



in + 1(1+1) J + ') 

and from below by the same expression with replaced by — ■ Using the 

identity 

2j + l 1 



1 



+ + J + l + l 1-3 



we find 



l-l 21-1 1-2 

E + E -2 E 

m=2 m=l+2 



1-2 



m=l 



m 



+ E ^( 



m 



-3 In/ + 



m=l 



n—l 



n+l+1 



lim 

n— +00 



lim In 

n— >oo 



E+E-2E £:+E^( 



_m=l m=2l+2 m=l+l 

{n-l){{n + l + l) 



m 



m=l+l 

+ ln(Z + l) + = In/ + 



Summing the expressions we get the upper bound in ([4.22|) ; the lower one is obtained 
in the same way. | 

Lemma 4.4 For any I large enough the interval (/(/ — I), /(/+!)) contains a point 
Hi such that ^{^/JTi) = . The number ni has the following properties: 

(i) l{l + l)-fii = 2l{\nl)-^{l + 0{l)), 

< In/ + C(l) for E<fii, 



m 



l{l+l)-E 
21-1 



0{1) for E>fii, 



l{l-l)-E 

w) finally, = -1 + 0((ln/)-^: 



Proof: Fix /. We have Z{kf-g{kf = (^Z{k) + {-iyg{k)) (z(A;) - (-l)'^(A:)) , and 
by the preceding lemmas these expressions equal 



2/±l 



27r/(/±l)-^ 2tt 



Inp + 0{1) 



The term with the minus sign is negative in (/(/ — 1), /(/ + 1)) provided / is large 
enough. The other term is sign changing for large / so it has has a root. In view of 
Lemma 4.1 there is just one fii such that 



2/ + 1 



In/ + = ; 



(4.23) 
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this proves (i). The relation ( [4 .231) yields also the next two claims: 

2/ + 1 . 2/ + 1 



and 



KM. 

21-1 
l{l-l)-E 



< 



E - l{l + 

2/-1 2/-1 



In/ + C(l) 



< 



f^i- 1(1-1) 2/(l-iigil)) 



Finally, (iv) follows from (ii) and (in) which yield 



i^(ln/ + 0(l))-i-ln/ + 0(l) 



An 



(9(1). 



In/ + C(l) 



(4.24) 



(4.25) 



and 



^?(VJ^) = ^in/ +0(1) 



Now we are in position to prove the main result of this subsection: 

Theorem 4.5 Let /sT^ := ^ \ (/";-^(0(lnO"^/^«+^(0(lnO"^) ; where e{-) is 
a positive strictly increasing function which tends to oo and obeys the inequality 
\e{x)\ < xliax for x > 1 . Then there is a positive c such that the transmission 
probability satisfies the bound 



\t{k)\^ < ce{iy 



(4.26) 



for A;^ G -R'e n (/(/ — I), /(/ + !)) and any I large enough. On the other hand, there 
are sharp resonance peaks localized on , 



\t{^{)\''= l + 0((ln/)"^ 



(4.27) 



as / ^ oo 



Proof: In the first part we are going to estimate the modulus of the numerator in 
(|4.2(]|) from below. The leading term at high energies is the one with A(fc) ; we shall 
estimate it on the set K^,] i.e., away of the zeros of the coefficient. Consider the 
neighborhood 



// := 



/ii-2/(ln/)"\/ii + 2/(ln/)"^ 



of /iz on which the following estimates are valid: 



^ {Z\k) + {-l)'g\k) 

1-2 oo 

E+ E 

m=0 m=l+l 

2(2/ + l) 



E 



2m + l 



„=o rn{m + l)-E 
2m+l 



l+m 



> 



> 



(m(m+l) -E)2' 
2(2/ + l 



1+ -1 



\l+m\ 



+ 



2(2/ + l) 



(/(/ + !)- E)^ 



(1(1 + 1) -E)^ - (/(/ + !)- /x, + 2/(ln/)-i) 

2(2/ + l) ^ (\nlf 

(2/(ln/)"i(2 + C(l))' ~ I 
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for some Ci > . Combining this result with Lemma 4.1 we are able to say how are 
the factors constituting A{k) separated from zero for E e (/(/ — 1), /(/ + 1)) which 
does not belong to the interval — e(/)(ln/)~^, {fii+e{l)(\ial)~'^) C Ii : 



Z{k) + {-iyg{k) > e{l){lnl)'\il-\\nlf = cirhil) . 



(4.28) 



Now we begin estimating the modulus of the transmission amplitude from below. 
In the following, Cj always means a positive constant. Using the expression ( [4.20| ) 
we get a simple lower bound, 



m\ < 



4:7ip 



l-27rZ(fc) 
k\9{k)\ 



(4.29) 



is obtained by neglecting the imaginary part of the denominator. First we shall 
show that 

^ - -(^) <«°) 

as / oo . Using Lemmas 4.2 and 4.3 we estimate this expression by 

21+1 I 21-1 



ln/ + 0(l) 
kK 



1 



2nk 



l(l+l)-E l{l-l)-E 



21+1 



21-1 



l{l+l)-E 



l{l-l)-E 



< C2. 



where in the second term we have used the explicit lower bound on \g{k) \ , estimated 
the modulus of the fraction by ^ < 5 , and employed finally k = a/E = / + 0{1) . 
For the second term in the denominator of ( [4.29| ) we shall show that 



'A(A;) 



{k-^-Akp^) > c^eil) 



holds for all Eg {l{l-l),l{l+l))\{pi-e{l){\nl)-'^ , pi+e{l){\nl) 



-2^ 



(4.31) 
which will give 



ill)with c:= (|^)' 
Lemmas 4.4(ii) yield 



TT 



k-^-Akp^ 

> ATT^p'^k 

e{l) 



. Consider first the case E < pi , when ( [4.28| ) together with 
Z{k) + {-iyg{k) Z{k)-{-iyg{k) 



> c^k ■ 



I 



J_ ( 21+1 
4n V 1(1+1)- 

1- (2pfc)-2 
> C5e{l) , 



21-1 



-E 



l{l-l)-E 



cii-Hi)l 


21-1 


- ln/ + 0(l 




l(l-l)-E 


) 




1 21- 
47r 1(1-1) 


1 1 
-E 47r 


ln/ + C(l) 





where in the last step we used again k = I + 0{\) . Let further E > pi . We divide 
the argument into two parts. First we suppose 



2/ + 1 
l{l + l)-E 



< 2 In / ; 



29 



then 



71 



Z{k) + {-^y9{k) Z{k)~{-iyg{k) 



21+1 



J_ 

47r V l{l+l)~E 



21-1 



l{l-l)-E 



+ Oil) 



cil-^ 


s{l) 


ln/ + C(l) 1 




ln/ + C(l) 





> An^p^k 1- {2pk) 

by ( [4 .281) and Lemma 4.4(iii). On the other hand, if 

2/ + 1 



> Cecil) 



lil + l)-E 

the same expression is bounded from below by 



> 2 In/ 



(4.32) 



Airy^k 1- i2pk) 



-2 



1 21+1 : Q/^^ 

1(1+1)- E ^ ^y^j 


^ In/ + 0(1) 




1 21+1 , 
47r l{l+l)-E ^ ^y^^ 





> cjk 1- (2pA;)-2 In/ > C8e(/) , 

where the denominator and the first term in the numerator have been estimated by 
means of ( |4.23| ) and ( [4.32|) , in the second term we have neglected one of the two 
terms of the same sign, and the last inequality follows from the fact that xlnx > 
eix) by assumption. To conclude the proof of (^4.26|) , it is sufficient to put C3 := 
min{c5, C6, cs} in (|4.31|) . 

The rest is easier; the existence of resonance peaks at which the sphere is almost 
transparent follows from the relations 



\t{f^i)\ 



1 



i + o((in/)-^; 



zi^{) 
) ai^M) 

(l + C»((ln/)-^ 



i + 



1 



1 



-1 



(ln/ + 0(l)) 



-1 



-1 



l + 0((ln/)-i) 



where we have used Lemma 4.4(iv). | 



4.4 A "bubble" array 

Having proved this we illustrate the behavior by numerical results. We find inter- 
esting the behavior of |t(A;)p. The curve oscillates almost regularly, but with the 
amplitude spaning from small values to unity for sufficiently large k. The lower 
enveloping curve of |t(fc)P behaves as iEliiE)'^ as we proved above. However, the 
behavior of |t(/c)P resembles that of ^'-interaction, if an smoothed (locally averaged) 
curve is compared, as already mentioned above. This suspicion can be supported 
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by Fig. 7; it compares an averaged of a "bubble" on the line and of 

(^'-interaction. The averaging is done over ten neighboring peaks of a given point. 
The strength of the ^'-interaction is chosen so as to reach the same value of 
at a distant k. These two curves seem to have the same "asymptotic" behavior. 

Fig. 8 shows the band spectrum of an infinite array of bubbles on the line. What 
is different from analogical spectra of loop or comb arrays is the concentration of 
bands to small values p (note the logarithmic scale). 



A A direct derivation of the comb— graph S— matrix 

To illustrate the advantages of the formulas ( p.31| ) and ( p.32|) , we present two other 
ways to derive the comb-shaped graph S-matrix. 

The first proof: By brute force, looking directly for the generalized eigenvector 
of if TV at the energy k"^ . One takes for the line component of the wavefunction the 
following Ansatz, 



^ikx _^ ikx 



-ikx 



ikx 



X <0 

L(j-l) <x<Lj 
X > L{N-1) 



(A.l) 



for j = 1, . . . , N—1 , while on the appendices we have multiples of the above specified 
solution to ( p.6|) , Us = Ps^i , s = 1, . . . , N . We have to find the coefficients for which 
such a vector belongs locally to the domain of the Hamiltonian. 

Substituting from ([A.l| ) to (|3^) we get a system of linear equations; denoting 



Co := 1, /„ := 0, 
we can write it concisely as 



Cat 



^N, fo 



rN, 



(A.2) 



-J 



= 0, (A.3) 

e,+,be^ + f,+,be-^ + (3,+,{cu'M - = 0, (A.4) 

{ik-d)e^ej+i + {ik + d)e'^ fj+i + ike^ej + ike~^ fj + (3j+ibUf,{0) = (A.5) 



for j = 0, . . . , N —1 . First we express ej from ( |A.3| ): by induction we check the 
relation 



(e 



(A.6) 



Using ( |A.4]) we exclude ej+i from ( |A.5| ), which yields 

{2ik-d)e~^fj - de^e^ - 2ike^^ fj+i + Pj+ibu'f,{0) 



0. 



(A.7) 
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Next we exclude Pj from (|A.4|) and (|A.7|) ; substituting then for Cj from (|A.6|) we 

get 



G 



-fiG + fo[G-F] = F, 
s-^f^-e^Fil-e')\J2fi^^A-^'fo = 0, (A.8) 



where j = 1, . . . , N—1 and F, G are the quantities defined above. The last system 
can be transformed in such a way that its matrix is tridiagonal; to this end we add 
the j-th equation to the (j — l)-th one multiplied by e . The resulting system has 
the following form, 



~hd{cu[-ue){0) + MG~F] 
-Gf,+, + [(1 + e')G + is'-l)F]f, - Gs'f,^, 



0, 



(A.9) 



where j = 2, . . . , N — 1 . In terms of (3, z this can be rewritten as 



/o(l + 2/5)-/i = -ifi, 

£~^fj+i - '^zfj + efj_i = , J = 1, 
-2zfN-i + efN-2 = ; 



(A.IO) 



we have taken into account here that /a? = . Now it is straightforward to find 
fo = as a ratio of the corresponding determinants, = -Oj^V-^lv'' 5 with 



D], = -iPDm^i , Dl = {l+t(3)DN-i+eD 



N-2 , 



where 



D 



N 



-2z e-^ ... 

e -2z e-^ 















-2z e-^ 
e -2z 



and Dq := 0. The above determinants satisfy the recursive relation Dnj^2 + 
2zDn^i + D]\r = , so by induction we get Dj^ = {—l)'^U]\f{z) , i.e., the first 
formula in (p.8|) . 

In a similar way we find the transmission amplitude. We start again from ( |A.3D 
from which we express fj , 



N-l 



£ [ei+i- 



j=0,...,N-l 



Transforming the system we arrive at an analogy of ( A.IO ) in the variables Cj , 

—2zei + 662 = , 

eej+2-'2zefj+i + e~^ej = 0, j = 1, . . . , iV-2 , 
eAr(l + iP) - Cn-I = , 
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which has to be solved for = ■ A disadvantage of this argument is that it gives 
a httle insight into the way in which graph appendices contribute to the S-matrix. 

The second proof: By S-matrix factorization described in Section 2.2. The on- 
shell S-matrix of a single tooth attached at x = is 



S = \ ^ ^-2ik^^ I 5 (A.ll) 



where t := ti and r := ri are the corresponding amplitudes of a tooth attached at 
X = which can be adopted from The matrix ( |A.11| ) is unitary since t/r is 
purely imaginary. The mentioned factorization technique yields for the amplitudes 
tn and the following recursive relations 

e^^rt\ ttN /A io\ 

rN+i = rN + 5 , tN+i = 5 , (A. 12) 

which are equivalent to 

tirN+i-r^) = e'^'^rtN+itN , {l-€\rN)tN+i = ttjv • (A. 13) 



Substituting from (|3^ ) into (|A.13| ) and employing the relations between Chebyshev 



polynomials we check that and satisfy indeed (|A.12| ). Certain disadvantage 
of this argument is that the closed form of the amplitudes ( |3.8| ) has to be guessed 
or derived by other means. 
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Figure captions 

Figures 2, 5, 8 available upon request from tater@ujf.cas.cz 
Figure 1 A loop-graph scatterer in a magnetic field 

Figure 2 The relation of transmission probabilities of a finite array of loops 
{N — 6) and the spectrum of the corresponding infinite system, Li — 0.5, 
L2 = 1.5, e^y/2,ai^ -1, a2 = -2. 

Figure 3 A comb-shaped graph 

Figure 4 The relation of transmission probabilities of a single appendix (dotted 
line) , a finite array of appendices {N — 7) (full line) , and the spectrum of the 
corresponding infinite system (thick lines above). The parameters are b — 1.2, 
c = 0.5, d = 0, L = 1, and£= 1. 

Figure 5 The band spectra of infinite arrays of appendices. The upper figure is 
for b — 1, c — d — 0, and the lower one for b = 3, c = 0.5, and d = —11. The 
spacing is chosen ^ = 1 in both cases. The full lines show the role of perfect 
transmission (the lines in the bands) and total reflection (the lines in the gaps) 
for a single appendix. 

Figure 6 This figure shows positions of resonances for = 5 appendices and 
b = 0.4, c — d = 0, and L = 1.2. The dotted line separates the two sets 
of resonances. The resonances closer to the real axis have their origin in A^- 
fold eigenvalue for b — 0; increasing hfts the degeneracy and pushes the 
resonances from the real line. The rest of the resonances (below the dotted 
line) has the origin in the serial arrangement of individual scatterers. 

Figure 7 The transition probability of a sphere of unit radius on a line (p = 
0.01). The lower graph compares the asymptotic behavior of |t(/i;)P for 6'- 
interaction with suitably chosen strength and that of averaged transmission 
probability plotted in the upper graph. The averaging is done over ten neigh- 
boring peaks of each point. 

Figure 8 Band spectrum of an infinite "bubble" array. The spheres are of unit 
radius, the spacing is £ = 1 (upper figure) and i — 0.01 (lower figure), p is the 
contact radius. 
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